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Abstract 

This is the completion of an exploratory study of Compact lattice Quantum Elec- 
trodynamics with a weak four-fermi interaction and four species of massless fermions. 
In this formulation of Quantum Electrodynamics massless fermions can be sim- 
ulated directly and Finite Size Scaling analyses can be performed at the theory's 
chiral symmetry breaking critical point. High statistics simulations on lattices 
ranging from 8 4 to 24 4 yield the equation of state, critical indices, scaling func- 
tions and cumulants. The measurements are well fit with the orthodox hypothesis 
that the theory is logarithmically trivial and its continuum limit suffers from Lan- 
dau's zero charge problem. 
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1 Introduction 



This is the final paper in a series of simulation studies searching for a nontrivial 
ultra-violet fixed point in abelian gauge theories. Previous studies have provided 
simulation evidence that scalar QED and noncompact lattice QED with four 
species of massless fermions Q are logarithmically trivial. 

This work employs the xQED formulation of the model in which a weak four- 
fermi interaction is added to the standard action. This affords us two advantages 
over standard methods: 1. we can simulate massless fermions directly on the 
lattice and see how massless fermion charge screening affects the dynamics, and 
2. the four-fermi interaction separates the chiral transition of the model from 
its confinement/deconfinement transition which is a first order transition and is 
controlled by monopole condensation. 

The xQED action contains two couplings, j3 = 1/e 2 where e is the usual elec- 
trodynamic charge, and the four-fermi coupling, G = 1/A. The model's phase 
diagram contains two separate lines of transitions, one describing monopole con- 
densation which is first order and a line of second order chiral transitions. We 
will show that the line of second order chiral transitions describes a logarithmi- 
cally trivial continuum limit, presumably identical to the continuum limit of the 
noncompact xQED lattice model studied earlier [0. 

This will be accomplished by exploiting the fact that Finite Size Scaling (FSS) 
applies simply to xQED because it is formulated without a bare fermion mass. In 
the conventional lattice action the bare fermion mass explicitly breaks chiral sym- 
metry and introduces scale breaking. A nonzero fermion mass is needed so that the 
standard algorithm converges. By contrast, xQED allows us to study the physics 
of the critical point by doing simulations in its immediate vicinity and extracting 
information from FSS arguments without ever needing uncontrolled extrapola- 
tions to the chiral limit. One of the purposes of this paper is the illustration of FSS 
methods for lattice gauge theories with massless fermions. This represents new 
territory for lattice gauge theory. 

Of course, this work is not without its disappointments. We cannot simulate 
the model in a range of parameters where monopoles are relevant degrees of free- 
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dom in the theory's continuum limit. For this action, the line of monopole con- 
densation transitions is distinct from that of the chiral transitions and monopole 
condensation appears to be first order. It has been argued elsewhere that a sec- 
ond order transition where there is both chiral symmetry breaking and monopole 
condensation would be a natural scenario for a nontrivial form of continuum QED 
because screening effects coming from fermions could balance anti- screening ef- 
fects coming from monopoles. A future formulation of lattice monopoles and 
fermions will have to be developed to see if this possibility could work out. 

Since this paper is a continuation in a series we refer the reader to reference 
for additional background and formulas. Here our emphasis is on the new 
simulations, FSS methods and results. 

The paper begins by laying out the lattice theory's two dimensional phase di- 
agram and presenting some simulations which pin down its qualitative features. 
Then some highly accurate 16 4 and 24 4 simulations in the broken symmetry phase 
are presented and analyzed. This is the approach used in past studies which 
showed that the noncompact model is logarithmically trivial Q. The compact 
model should also be logarithmically trivial because the chiral transition occurs 
in a region of the model's phase diagram where monopoles are not critical. Un- 
der these conditions the differences between the compact and noncompact models 
should become irrelevant and the models should have the same continuum limit 
0. This conventional idea was not well supported by our first simulation study 
of the model which motivated us to do a more thorough job which we report 
here. The conclusions presented here are based on over ten times the statistics of 
earlier work. High statistics allowed the simulation program to reach thermody- 
namic equilibrium and produce more accurate ensembles of configurations which 
reduced the error bars substantially. The 16 4 and 24 4 data in the broken sym- 
metry phase can be fit with both power laws whose critical indices deviate from 
mean field theory, as reported in 0, or log-improved mean field theory. The log- 
improved mean field theory fits are somewhat better than the power law fits, but 
a conclusive result eludes us if we just use this subset of our simulations. How- 
ever, the combination of the broken phase results together with new simulations, 
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analyzed with FSS methods on lattices ranging from 8 4 to 16 4 , in the immediate 
vicinity of the model's critical point lead to stronger results: they favor the log- 
arithmically trivial fermionic field theory scenario. This is the new conclusion 
which will be presented below. 

This paper is organized as follows: In the next section we briefly review the 
formulation of the lattice action and define the parameters we use in our fits and 
simulations. In the third section we sketch the phase diagram and present a few 
simulations which led to it. Then we examine several points in the phase dia- 
gram along the line of chiral transitions and show that the data for the chiral order 
parameter in the broken symmetry phase is compatible with log-improved mean 
field theory. Next we review the relevant features of FSS and present several sec- 
tions of analysis of high statistics data sets using lattices ranging from 8 4 through 
16 4 . This is the most decisive analysis in this study. Its success depended crucially 
on the use of the xQED action and statistically large data sets which are several 
orders of magnitude larger than those used in typical lattice simulations of QCD, 
for example. FSS simulations at and near the critical coupling were essential to 
establishing the logarithmic triviality of this model. 

2 Formulation 

The lattice Action of compact xQED, where the gauge symmetry is interpreted as 
a compact local U(l) symmetry, following Wilson's original proposal |6|, reads 



S = Hx)(M xy +D xy )^y)+— * 2 (*) + ^ E (l-cos(F^(x))), (1) 



F^(x) = e il (x) + 6 u (x + fi)+6_^x + fi + u) + 6_ u (x + v), (2) 



X 



where 




(3) 
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The auxiliary scalar field a is defined on the sites of the dual lattice x Q, and the 
symbol (x, x) denotes the set of the 16 lattice sites surrounding the direct site x. 
The factors e ±t9 ^ are the gauge connections and rj^x) are the staggered phases, 
the lattice analogs of the Dirac matrices. ^ is a staggered fermion field and m is 
the bare fermion mass, which will be set to zero. Note that the lattice expression 
for F^v is the circulation of the lattice field 6^ around a closed plaquette, the gauge 
field couples to the fermion field through compact phase factors to guarantee local 
gauge invariance and cos enters the action to make it compact. 

It will often prove convenient to parametrize results with the inverse of the 
four-fermi coupling, A = 1 /G, and the inverse of the square of the gauge coupling, 
P = 1/e 2 . 

The global Z 2 chiral symmetry of the Action reads: 



tj}{x) -> {-1) x1+x2+x3+x4 i/j{x) (5) 
ij(x) -> -^( x )(-l) xl+x2+x3+x4 (6) 
a -> -a. (7) 

where (^—i Y 1 + x2 + x3 + xi [ s the lattice representation of 75. 

Interesting limiting cases of the above Action are: (i) the Z 2 Nambu— Jona- 
Lasinio model with no gauge fields, set e 2 to zero here, which has a logarithmi- 
cally trivial chiral phase transition at nonzero G; (ii) the compact QED model with 
no four-fermi interactions, whose first order chiral phase transition is coincident 
with its first order monopole condensation transition near j3 = 1/e 2 ~ 0.89(1) 
for four flavors [9|; and (iii) the G — > 00 limit in which the fermions obtain a 
dynamical mass comparable to the reciprocal of the lattice spacing and therefore 
decouple, leaving quenched compact QED which has a first order transition at 
(3 = 1.011124(1) ifTOll . 

We refer the reader to earlier papers in this series Q for more motivation and 
details. We only emphasize new simulations and analyses here. 
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Figure 1: Phase Diagram of Gauged Compact U(l) Nambu Jona-Lasinio Model 

3 Overview of the Phase Diagram. 

The phase diagram for the model |3 1 is shown in Fig. [lj The monopole concentra- 
tion M and the chiral condensate, (^), label the three phases. Region I has chiral 
symmetry breaking in a condensate of monopoles, Region II has chiral symmetry 
breaking in a monopole free vacuum, and Region III is chirally symmetric in a 
monopole free vacuum. Past simulations suggest that the dashed line consists of 
second order transitions and the thick line consists of first order transitions Q. 

It is interesting to confirm that the dashed line of chiral symmetry breaking 
transitions in the upper reaches of Fig. dturns vertical and the four-fermi coupling 
alone breaks chiral symmetry at strong coupling, in agreement with [8|. This is 
shown in Fig. 2 where we confirm that the model breaks chiral symmetry for 
G = 1.0 no matter how small the gauge coupling e 2 = 1/(3. 

In reference Q we concentrated on the vertical line G = 1/1 A and showed 
that the line of first order monopole concentration transitions is clearly separate 
from the line of chiral symmetry transitions. We felt that it was important to 
simulate along another vertical line in the phase diagram to verify this result and 
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Figure 2: (a) vs. gauge coupling [3 = 1/e 2 at fixed four-fermi coupling G — 1.0. 



to see how sensitive the characteristics of the transitions are to the bare couplings. 
There are several questions we need to address. They include: (i) Does the order 
of the transition(s) change along the lines in the phase diagram? (ii) Do the critical 
indices change along the lines? There are related models where this is known to 
happen, such as in the quenched noncompact gauged Nambu-Jona Lasinio model 
ifTTl . It is unknown if such phenomena can occur in unquenched abelian models 
where fermion screening leads to the zero charge problem in perturbation theory 
and produces only perturbatively trivial models. 

In Fig. 3 we show the monopole concentration as the gauge coupling (3 passes 
through the phase boundary at G = 0.50 between regions I and II. The monopole 
concentration appears to have a discontinuous jump just below j3 = 0.935. 

Measurements of (a) in Fig. 4 also show the first order transition near j3 = 
0.935 and indicate that region II (see Fig. [TJ extends from this point to approxi- 
mately j3 ~ 0.96 along the vertical line at G = 0.50. 

The simulations reported in Q along the vertical line G = 1/1 A found a con- 
tinuous chiral transition at (5 — 1.393(1). It is interesting to check the consistency 
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Figure 3: Monopole concentration vs. gauge coupling (3 at fixed four-fermi cou- 
pling G = 0.50. 
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Figure 4: (a) vs. gauge coupling (3 at fixed four-fermi coupling G = 0.50. 
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Figure5: (a) vs. four-fermi coupling A = l/G at fixed gauge coupling j3 = 1.393. 



of this prediction with simulations approaching the critical point from a different 
direction in the phase diagram. Therefore, we simulated the model on 16 4 lat- 
tices along the horizontal direction, fixing (3 = 1.393 and varying the four-fermi 
coupling. The results are shown in Fig. 5. 

We see that Fig. 5, although not competitive in accuracy with the study at 
fixed four-fermi coupling and variable gauge coupling due to an apparently narrow 
scaling region, is indeed compatible with the earlier results: the critical point is 
again predicted to be at G = 1/1.4 and (3 = 1.393. The dashed curve in Fig.|5]is 
just meant to guide the eye to the horizontal line. 

4 High Statistics 16 4 Simulations along the Vertical 
Line A = 1.4 

The real focus of this series of simulations is to determine the quantitative nature 
of the line of chiral transitions in the phase diagram. We, therefore, developed 
a very fast, parallel version of the simulation code and made high statistics runs 




Figure 6: (a) vs (3 at four-fermi coupling A = 1.4. The fit has the critical index 

(3mac = 0.65(4). 



along the vertical line A = 1.4, complementing the exploratory runs reported 
earlier Q. In fact, between 10 x 10 6 and 20 x 10 6 sweeps of the Hybrid Molecular 
Dynamics algorithm were done at each of fifteen couplings j3 ranging from (3 = 
1.393 to j3 = 1.125. This is more than an order of magnitude greater statistics than 
those previously reported and represent between 100, 000 and 200, 000 trajectories 
of the Hybrid Molecular Dynamics algorithm lfi~2l at each coupling (the Monte 
Carlo time interval was chosen to be dt = 0.01 to keep systematic errors negligible 
in the molecular dynmics steps). The data on the broken chiral symmetry side of 
the transition is plotted in Fig.|6]and a power-law fit is included in the figure. (The 
data on the other side of the transition will be used in a FSS analysis below.) The 
power-law fit, (a) = A(j3 c - (3Y ma9 , is acceptable, x 2 /DOF = 1.7 with a critical 
point j3 c = 1.33(1) and a critical magnetization exponent (3 mag = 0.65(10). 

It is interesting that the critical index is larger than the mean field value of 
1/2, but the deviation from mean field theory is less than reported in our earlier, 
lower statistics, exploratory work Q. The high statistics of this run are making a 
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difference in the results. We recall from other works using U(l) gauge fields and 
the simplest Wilson action, that very long relaxation times are noted (1701. Our 
time correlation analysis suggests that the 20 million sweeps used here suffice and 
the error bars in the figure conservatively account for the correlations in Monte 
Carlo time. It is interesting that simulations of the noncompact model [2 | did not 
require such enormous statistics to achieve equally accurate results. 

Following our analysis of the noncompact model, however, it is interesting to 
attempt to fit the same data with the hypothesis of logarithmic triviality. As dis- 
cussed in previous work, including reference [ 81, the logarithms of triviality effect 
the scaling laws and equation of state differently for fermionic theories than for 
bosonic theories. In particular, in references EUHl, the leading order equation of 
state had the form (3 C — (3 = A(a) 2 (hi{l/ (a)) + B). Fits of this form accomodate 
the data of Fig. |6]very well, in fact, just like in the noncompact model. Fitting 
routines predict the parameter B = 0.84(8), A = 7.3(9) and (3 C = 1.33(1) with 
X 2 /DOF = 0.66. To emphasize the need for the logarithm here, we plot the quan- 
tity (P c — (3)/ {a) 2 against ln(l/(cr)) in Fig.|7J and show the fit. This makes the 
point that the simulation results are compatible with logarithmic triviality. In fact, 
they appear to rule out bosonic triviality fits which would have the logarithm of 
triviality in the denominator of the equation, like (3 c —f3 = A(a) 2 / (ln(l/ (a))+B), 
rather than the numerator. Our work, therefore, supports the analytic predictions 
of A. Kocic lfT3l . Perhaps, the fact that the x 2 /DOF from the log-improved 
mean field relation (which is a three-parameter function) is less than half the 
X 2 /DOF we get from the fit to the standard power-law relation (which is also 
a three-parameter function) is evidence that the data favor the triviality scenario 
over the interacting field theory scenario. More compelling and straight-forward 
evidence for triviality will be presented in the sections on FSS below. 

This result is not above criticism, however. As is clear from the figures, we 
are not able to simulate the model very close to the critical point on this lattice 
size without meeting uncontrollable finite size effects. For example, simulations 
closer to (3 C display tunneling between the Z 2 vacua through the unbroken phase 
and make reliable measurements of the order parameter impossible. We will turn 
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Figure 7: (/3 C — 0)/ (a) 2 vs ln(l/ (cr)) at Four Fermi coupling A = 1.4. 

to FSS simulations below to remedy this limitation and obtain a better estimate 

for (3 C . 

5 24 4 Simulations along the Vertical Line A = 2.0 

Next consider simulations on larger lattices in the broken symmetry phase. In 
Fig. [8] we show the raw data from simulations on 24 4 lattices at a somewhat weaker 
four-fermi coupling G = 0.5. These larger lattices allowed us to run simulations 
closer to the critical point without suffering from large finite size effects. The 
algorithm was also somewhat better behaved at weaker four-fermi coupling G = 
0.5. Of course we were not able to amass as high statistics in this case: one million 
sweeps per coupling were accumulated. 

The raw data in Fig. [U is fit with a simple power, (cr) = A{[3 C — (3Y maa , and the 
parameters A = 2.8(9), (3 C = 0.952(1), (3 mag = 0.77(10), were determined with 
X 2 /DOF = 1.18. We again see that the best power law predicts a critical index 
Pmag higher than that of pure mean field theory, as in reference Q. However, 
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Figure 8: (a) vs (3 at Four Fermi coupling A = 2.0. The fit has the critical index 

Pmag = 0.77(10). 

the data is also well fit with the hypothesis of triviality. Fits to the form (3 C — 
(3 = A(a) 2 (ln(l/(a)) + B) gave the parameters A = 0.21(10), (3 C = 0.951(1) 
and B = -0.061(1) with a fine quality of fit \ 2 /DOF = 0.78. Following our 
presentation above, we plot {(3 C — (3) /(a) 2 against ln(l/(er)) inFig.|9j which shows 
the importance of the logarithms and the deviation from pure mean field theory. 
As in Sec. 4 we find that the log-improved triviality fit is preferred to the power 
law form. 



6 Results near the critical coupling 

6.1 Background on Finite Size Scaling 

In order to study the critical behavior of our theory arbitrarily close to the critical 
coupling we used FSS methods. FSS techniques first developed by Fisher lfT4l 
are important tools used in the determinations of critical exponents near second 
order phase transitions. The critical behavior of a system in the thermodynamic 
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Figure 9: (f3 c — 0)/ (a) 2 vs ln(l/ (a)) at Four Fermi coupling A = 2.0. 

limit may be extracted from the properties of finite size systems by examining the 
size dependence of the singular part of the free energy density. According to FSS 
theory, for dimensionality d less than the upper critical dimension d c , the singular 
part of the free energy is described phenomenologically by a universal scaling 
form, 

F s (t, m, L) = L~ d F(tL 1/u , mL ( ^^ )/v ), (8) 

where m is the fermion bare mass and t = (/3 c —j3). The critical exponents v, f3 mag 
and 7 are all the thermodynamic values for the infinite system. Scaling forms for 
various thermodynamic quantities can be obtained from appropriate derivatives of 
the free energy density. On a finite volume and with the fermion bare mass set to 
zero, the direction of symmetry breaking changes over the course of the run so the 
chiral condensate averages to zero over the ensemble. Another option is to intro- 
duce an effective order parameter E = (\cr\), which in the thermodynamic limit 
is equal to the true order parameter (a). The FSS scaling form for E determined 
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from Eq.[5]is 

£ = L- Pmaalu f u {tL l,u ). (9) 

A standard method to measure (3 C for a second order transition is to com- 
pute the Binder cumulant [fT31l for various system sizes. On sufficiently large 
lattices where subleading corrections from the finite lattice size L are negligible, 
the Binder cumulant Ub((3, L), defined by 

U^l- 1 -^, (10) 

is given by U B — fsLitL 1 ^) and, therefore, at f3 c it becomes independent of L. 

Another quantity of interest is the susceptibility x which, in the static limit of 
the fluctuation-dissipation theorem, is 

X = lim V[(a 2 ) - (a) 2 }, (11) 

where V is the lattice volume. For finite systems this expression leads to the 
following finite-lattice estimates for x- 

Xi = V(a 2 ) f3>p c , (12) 

Xc = V[(a 2 )-(\a\) 2 ] (3<f3 c , (13) 

where the subscript c stands for "connected." Both relations should scale at criti- 
cality like L 1 ^ . Furthermore, the maxima of Xc m the scaling region should also 
obey xr ak ~ Lljv ■ 

Furthermore, logarithmic derivatives of (|cr| n ) can give additional estimates 
for v. It is easy to show that 

'(W\ n P) 



D n = ^-ln(\a\ n ) 



-<P) 



(14) 



where P is the plaquette, has a scaling relation 

D n = L^f Dn (tL^)- (15) 
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Other related quantities useful in determining the critical exponent v can be 
defined from logarithms of derivatives of (a n ) |[T6l . In our analysis we will con- 
sider 

Q^2[a 2 ]-[a% (16) 

where 

d(a n ) 

One can easily show that 

Q~-hiL+Q{tL l i v ). (18) 

The above FSS relations rely on the traditional FSS hypothesis that in the 
vicinity of the critical coupling the behavior of the system is determined by the 
scaled variable L/£, where £ is the correlation length. The standard FSS hypoth- 
esis fails for d > d c . A modified hypothesis for O(N) -symmetric $ 4 theories in 
four dimensions was proposed in ifTTll . where it was shown that in the vicinity of 
the critical coupling the actual length of the finite size system is replaced by its 
correlation length £l(0) oc L(lnL)*, independent of N. However, as shown in 
1 8 ] and demonstrated numerically in earlier sections of this paper, the logarithmic 
triviality in fermionic field theories such as the NJL model and QED is manifested 
in a different way from the triviality in purely bosonic theories. The logarithmic 
corrections in the scaling relations of QED are expected to be in the denomina- 
tor of the scaling functions, whereas in $ 4 theory they are in the numerator. The 
same is expected in the FSS relations of various thermodynamic quantities. By 
generalizing the Privman-Fisher ansatz for the scaling relation of the singular part 
of the free energy, the scaling function for F s (t, m , L) becomes 

F„(t, mo, L) = L~ d T{tL 2 \^ v L, m L 3 \n q L). (19) 

Consequently, the log-improved FSS relations for the effective order parameter 
and the susceptibility obtained from appropriate derivatives at zero fermion mass 
are, 

£ = L'^Lf^tL^L), (20) 
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X = L 2 \^Lf x {tL 2 WL). 



(21) 



In the next sections we will present our attempts to extract the critical expo- 
nents and to check consistency with log-improved scaling. We studied the FSS 
behavior of several observables in order to compare different results and reach 
conclusions. 

6.2 Analysis of the Binder cumulant 

In the vicinity of (3 C we can expand the Binder cumulant and find 



By fitting this relation to data in the range 1.2875 < (3 < 1.375 and lattice sizes 
L = 8,. .,16 we found v = 0.49(10), (3 C = 1.356(7) and U* = 0.153(26) with 
X 2 /DOF = 1.6. The measured value of the critical exponent v is consistent with 
the mean field prediction v = 1/2, although the error is relatively large. More 
measurements of v with better precision are presented below. The location of the 
critical point [3 C = 1.356(7) refines the estimates found from the broken symmetry 
fits given earlier. Since those studies were done far from f3 c , a small discrepancy 
is not surprising. 

6.3 Analysis of effective order parameter 

In this section we discuss the FSS analysis for the effective order parameter S. We 
fit our data to E 0C = aLr Pma 'l v for three values of f3 = 1.325, 1.350 and 1.393, 
which are close to the value (3 C = 1.356 extracted from the analysis of Ub(L, (3). 
The results are presented in Tableland plotted in Fig. [TTJ We conclude that (3 C 
is close to 1.350 and the ratio f3 mag /v = 1.14(6) at this value of (3 is also close to 
the mean field result (3 mag / v = 1. 

We also fit all the data in the vicinity of the transition to a single scaling func- 
tion obtained from the Taylor expansion of Eq. 9, up to a linear term, 



U B {(3 : L)^U, + U l {(3 c -f3)L 1 /\ 



(22) 



(23) 
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Figure 10: Binder cumulant versus coupling for different lattice sizes. 

After fixing (3 C = 1.356 we obtained the values of (3 m ag/v and v. When all the 
three values of the coupling {(3 = 1.325,1.350 and 1.393) for L = 8, ...,16, 
and an extra coupling (3 = 1.375 for L = 16, are included in the fit we get 
Pmag/v = 1.13(4) and v = 0.48(8) with x 2 /DOF = 1.4. For the same data 
set and for fixed (3 C = 1.330 (which is the value of the critical coupling obtained 
from the broken phase analysis) we get (3 mag /u = 0.95(5) and v = 0.56(8) with 
X 2 /DOF = 1.3. 

In order to check whether our results are consistent with log-improved mean 
field scaling, we fit £ to Eq.|20l The results are summarized in Tableland plotted 
in Fig. [T2| The best fit is at (3 = 1.350 with p = —0.34(14). The negative sign of 
p is consistent with the scenario of log-triviality in fermionic field theories [13]. 

6.4 Analysis of susceptibility 

First, we fit \i (Eq.fT2l as a function of L at different values of (3 = 1.325, 1.350 
and 1.393. The results are displayed in Fig.[l3]and Table|3] It is clear from these 
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Figure 11: S vs. L for different values of (3. 



results that the critical coupling is close to (3 = 1.350 in agreement with analyses 
presented in previous paragraphs. The measured value ofj/v = 1.76(14) is close 
to its mean field value of 2. We also fit xi m the vicinity of the transition to the 
linear expansion of x{L, (3) around (3 C . 



Xi(L,f3)^{c 1 + c 2 (p c -P)L^}L^. 



(24) 



After fixing the critical coupling to the value extracted from the broken phase 
analysis (J3 C = 1.330) we get = 0.66(9) and 7 /z/ = 2.15(9) with x/DOF = 1.3, 
whereas after fixing f3 c to the value extracted from the U b analysis we get u = 
0.51(7) and r y/v= 1.87(8) with the same x 2 /DOF as before. Furthemore, we 
fit the data to xi — aL 2 (\n 2p L) in order to check whether the data are consistent 
with log-improved mean field scaling. The results displayed in Table |5] provide 
good evidence that the log-improved scaling relation describes the data well and 
that (3 C ~ 1.350. The measured value p = —0.28(16) is compatible with the result 
extracted from log-improved fits of E discussed in the previous section. 
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Figure 12: Fits to LE = a ln p L for different values of (3 near the transition. 

We repeated the above analysis for the connected susceptibility \c- The re- 
sults presented in Table |5] indicate clearly that j3 c is close to (3 ~ 1.350 with the 
value 7/1/ = 1.70(11) close to the mean field result. The results of fits to the log- 
improved FSS relation are summarized in Table |6] Again the results provide sig- 
nificant evidence that (3 C is close to 1.350 and p = —0.36(13) which is consistent 
with previous measurements of p. We also fitted the peaks of Xc to Xc° ak ~ 
and got 7/v = 1.99(16) in good agreement with the mean field prediction. The 
data and the fitting function are shown in Fig. [141 

6.5 Analysis of Dj and Q 

To make a further check of our results for the exponent v, we studied the finite 
size scaling properties of logarithmic derivatives of (|cx| n ) defined in Eq.[l4j We 
fit Dj ~ L l l v for j = 1,2, and 3 at different values of the gauge coupling and 
for L = 8, 14. We note that the data generated on 16 4 lattices were noisy and 
therefore could not be included in the fits. The results for the exponent v and the 
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Figure 13: Susceptibility xi vs. L for different values of /3. 
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Figure 14: Peaks of connected susceptibility \c vs. L. 
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quality of each fit are shown in Tables 1718191 ithe values of v are in very good 
agreement with the mean field prediction v = 0.5. Although for these values 
of gauge coupling v do not have a significant dependence on j3, the qualities of 
the fits indicate that j3 c ~ 1.350 which is in agreement with results presented in 
previous paragraphs. Our attempts to fit the data to log-improved FSS mean field 
scaling laws Dj ~ L 2 ln p L did not give any signal for p. 




Figure 15: D u D 2 and D 3 vs. L at = 1.350. 

Furthermore, we fit all the data {(3 = 1.325, 1.350, 1.393) to the linear expan- 
sion of the FSS relation 

L) ~ Clj L l l» + c 2j {(3 c - /3)L 2 ' V , (25) 

with fixed j3 c = 1.356. The results are summarized in Table fTOl and are in good 
agreement with the mean field result although the quality of each fit is not as good 
as before. 

Finally, a fit to the linear expansion of the observable Q (defined in Eq.[T8l 

Q w - In L + d + c 2 (P c - (3)L 1/V (26) 
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for/3 = 1.300,..., 1.393, L= 8, 14 and fixed & = 1.356 gives v = 0.52(3) and 
X 2 /DOF = 1.5, in good agreement with all the measurements for v presented in 
previous paragraphs. 

7 Conclusions 

We have presented evidence for the logarithmic triviality of the chiral transition 
line in compact QED with four species of fermions. A weak four-fermi interaction 
was employed in the action so that massless fermions could be simulated directly, 
thus avoiding the need to extrapolate raw data to the chiral limit. xQED allowed us 
to use the simplest single variable finite size scaling fits to the data. Since the four- 
fermi interaction is irrelevant in four dimensions and since the full model xQED is 
found to be logarithmically trivial, this study constitutes strong evidence that the 
continuum limit of the standard compact lattice QED model is also logarithmically 
trivial. 

Finite size scaling proved to be an effective approach to deciding the physics 
issues inherent in these models. The FSS analyses provide strong evidence that 
the critical exponents are the mean field theory ones. Especially in the case of 
v the analyses of various observables show in a consistent manner that v is very 
close to 0.5. The broken phase data at different values of the four-fermi coupling 
also favor fermionic log-improved mean field scaling over the interacting field 
theory scenario. 
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Table 1: Results from fits to £ = aL~ Pma ^ v . 



p 


Pmag / ^ 


a 


X 2 /DOF 


1.325 


0.84(8) 


0.46(9) 


2.1 


1.350 


1.14(6) 


0.84(12) 


0.16 


1.393 


1.45(5) 


1.48(20) 


4.6 



Table 2: Results from fits to £ = aL' 1 \n p L. 



p 


a 


P 


X7DOF 


1.325 


0.49(8) 


0.37(19) 


2.2 


1.350 


0.80(10) 


-0.34(14) 


0.20 


1.393 


1.26(14) 


-1.06(13) 


5.4 



Table 3: Fits for X i = aL^' v '. 



p 


a 


l/v 


x 2 /dof 


1.325 


0.29(10) 


2.36(14) 


2.0 


1.350 


0.94(30) 


1.76(14) 


0.12 


1.393 


2.7(7) 


1.18(11) 


4.0 



Table 4: Results from fits to Xi = aL 2 {\n 2q L). 



p 


a 


P 


xVdof 


1.325 
1.350 
1.393 


0.33(10) 
0.86(23) 
2.06(45) 


0.41(17) 
-0.28(16) 
-0.97(13) 


2.2 
0.13 
4.7 


Table 5: Results from fits to Xc = aL 1 ^ . 


P 


a 


l/v 


X 2 /DOF 


1.325 
1.350 
1.393 


0.056(15) 
0.36(9) 
0.80(18) 


2.51(11) 
1.70(11) 
1.23(9) 


31.3 
0.9 
4.5 
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Table 6: Results from fits to Xc = aL 2 (\n 2p L). 



p 


a 


P 


xVdof 


1.325 


0.11(3) 


0.30(17) 


45.6 


1.350 


0.33(7) 


-0.36(13) 


0.9 


1.393 


0.62(12) 


-0.91(11) 


5.1 



Table 7: Results of fits to D 1 ~ L l l 



p 


V 


rVDOF 


1.325 


0.46(4) 


i.i 


1.350 


0.48(4) 


0.4 


1.393 


0.47(6) 


3.3 



Table 8: Results of fits to D 2 ~ L l l 



p 


V 


x 2 /dof 


1.325 


0.47(4) 


1.0 


1.350 


0.47(5) 


0.3 


1.393 


0.48(7) 


2.8 



Table 9: Results of fits to D z ~ L 1 / 



p 


V 


X 2 /DOF 


1.325 


0.49(4) 


1.0 


1.350 


0.47(5) 


0.3 


1.393 


0.51(7) 


2.3 



Table 10: Results of fits to eg. 1251 





V 


xVdof 


D l 


0.49(3) 


2.5 


D 2 


0.51(3) 


2.0 


D 3 


0.51(3) 


1.9 
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